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ABSTRACT 

We have begun a systematic numerical study of the nonlinear growth of nonaxisymmetric perturbations 
during the ambipolar diffusion-driven evolution of initially magnetically subcritical molecular clouds, 
with an eye on the formation of binaries, multiple stellar systems and small clusters. In this initial 
study, we focus on the m = 2 (or bar) mode, which is shown to be unstable during the dynamic collapse 
phase of cloud evolution after the central region has become magnetically supercritical. We find that, 
despite the presence of a strong magnetic field, the bar can grow fast enough that for a modest initial 
perturbation (at 5% level) a large aspect ratio is obtained during the isothermal phase of cloud collapse. 
The highly elongated bar is expected to fragment into small pieces during the subsequent adiabatic 
phase. Our calculations suggest that the strong magnetic fields observed in some star-forming clouds 
and envisioned in the standard picture of single star formation do not necessarily suppress bar growth 
and fragmentation; on the contrary, they may actually promote these processes, by allowing the clouds to 
have more than one (thermal) Jeans mass to begin with without collapsing promptly. Nonlinear growth 
of the bar mode in a direction perpendicular to the magnetic field, coupled with flattening along field 
lines, leads to the formation of supercritical cores that are triaxial in general. It removes a longstanding 
objection to the standard scenario of isolated star formation involving subcritical magnetic field and 
ambipolar diffusion based on the likely prolate shape inferred for dense cores. Continuted growth of the 
bar mode in already elongated starless cores, such as L1544, may lead to future binary and multiple star 
formation. 

Subject headings: binaries: formation — ISM: clouds — ISM: magnetic fields — MHD — stars: 
formation 



1. INTRODUCTION 

A basic framework has been developed for the forma- 
tion of low-mass stars in relative isolation (Shu, Adams, & 
Lizano 1987). It involves gradual condensation of dense 
cores from strongly magnetized, background molecular 
clouds through ambipolar diffusion, followed by dynamic 
core collapse to form stars. Quantitative studies based 
on this by now "standard" scenario have been carried 
out by many authors, with increasingly sophisticated in- 
put physics (Mouschovias & Ciolek 1999). However, with 
few exceptions (Indebetouw & Zweibel 2000; Boss 2000), 
axisymmetry has been adopted. The adopted symmetry 
precludes a detailed investigation of cloud fragmentation, 
generally thought to be a necessary step in the forma- 
tion of binary and multiple stellar systems; it is in such 
systems that most stars are found. We seek to improve 
this situation by removing the restriction of axisymmetry 
and systematically investigating the effects of the strong 
magnetic fields envisioned in the standard picture of single 
star formation on fragmentation and their implications on 
binary and multiple star formation. 

Fragmentation of nonmagnetic clouds has been studied 
extensively over the years, mostly through numerical sim- 
ulations (Bodenheimer et al. 2000, and references therein). 
In the canonical case of a spherical cloud, the fragmenta- 
tion is mainly controlled by the ratios of the cloud thermal 
and rotational energy to the gravitational energy, a and /3, 



although the distributions of density and angular momen- 
tum also play a role. Molecular line observations (Good- 
man et al. 1993) suggest that star-forming cores of molec- 
ular clouds are far from being rotationally supported, with 
a typical (low) value oi (5 ^ 0.02. For such slowly rotating 
cores, the criterion for fragmentation is roughly a <^ 0.5, 
if the cores are idealized as spheres of uniform density and 
rigid rotation (Tsuribe & Inutsuka 1999). More realis- 
tic, centrally peaked density distributions would lower the 
critical a for fragmentation. Such distributions are com- 
monly inferred in the well-observed starless cores such as 
L1544 (Andre, Ward-Thompson, & Barsony 2000; Evans 
et al. 2001). The central concentration, if taken as an ini- 
tial condition, tends to make fragmentation more difficult 
(Boss 1993). 

Fragmentation of strongly magnetized clouds with mass- 
to-flux ratios below the critical value (27r)^^G'~^/^ (i.e., 
subcritical) is expected to be quite different (Galli et al. 
2001). Through magnetic braking, the field controls the 
angular momentum evolution of the cloud, until a super- 
critical dense core forms (Basu & Mouschovias 1994); it 
thus determines the value of /3 of the core. More impor- 
tantly, the field can provide most of the cloud support 
against self-gravity and thus allow a cloud to have an arbi- 
trarily small thermal energy compared to the gravitational 
energy and still be in a mechanical equilibrium initially. 
The low value of a is a key ingredient of fragmentation. 
On the other hand, the presence of a strong magnetic field 
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tends to stifle cloud fragmentation, if the field and mat- 
ter are well coupled (Phillips 1986). Indeed, for a frozen- 
in subcritical field, fragmentation is suppressed altogether 
(Nakano 1988). Fragmentation can in principle resume in 
the part of cloud that has become magnetically supercriti- 
cal (Boss 2000), through ambipolar diffusion. It is however 
difficult to predict a priori how the firagmentation should 
proceed, if at all, in a cloud that is only partially coupled 
to a strong magnetic field. 

Langcr (1978) showed through linear analysis that, in 
a lightly ionized medium such as a molecular cloud, the 
Jeans instability is not completely inhibited no matter how 
strong the magnetic field is; rather, it grows on an ambipo- 
lar diffusion, rather than dynamic, timescale. Li (2001) 
followed the nonlinear evolution of a set of magnetically 
subcritical, super-Jeans mass clouds assuming axisymme- 
try and found that either a dense, supercritical core or 
ring forms as a result of ambipolar diffusion. The super- 
critical ring is expected to fragment readily into a number 
of dense cores, producing perhaps a small group of stars. 
The ring breakup will be studied elsewhere. In this Let- 
ter, we shall concentrate on the more subtle problem of the 
fragmentation of magnetized single cores, with an eye on 
forming binaries and multiple stellar systems. As is com- 
mon with numerical studies of cloud fragmentation and 
binary formation, we focus on the growth of m = 2 (or 
bar) mode (e.g., Nakamura & Hanawa 1997; Bodenheimer 
ct al. 2000). In § 2 we describe oiiv forniidation of the 
problem of nonaxisymmetric cloud evolution. Numerical 
results on bar formation are presented in § 3, and their 
implications on binary and multiple star formation as well 
as the observed elongation of dense cores are discussed in 
§4. 

2. PROBLEM FORMULATION 

An ordered, subcritical magnetic field allows cloud ma- 
terial to settle along field lines into a disk-like geometry. 
We assume that the disk is in hydrostatic equilibrium in 
the vertical direction, and adopt the thin-disk approxima- 
tion (Li 2001). The cloud evolution is followed in the x-y 
plane of a Cartesian coordinate system {x,y,z). Outside 
the disk, the magnetic field is assumed to be current free 
and uniform far from the cloud. An isothermal equation 
of state is adopted (unless noted otherwise), so that the 
pressure and surface density are related through P = c^E, 
where Cs is the (effective) isothermal sound speed. The 
governing equations are the nonaxisymmetric versions of 
equations (1) through (7) of Li (2001). 

The initial conditions for star formation are not well de- 
termined either observationally or theoretically. Following 
Basu & Mouschovias (1994) we prescribe an axisymmetric 
reference state for our model clouds with the distributions 
of mass and magnetic field given by 

Sref(a;,y) = ^"'"''^ and B^^^ei{x,y) = B^^^, 

[1 + r'^/r^] 

(1) 

where r = (x-^ + y'^Y^'^- The background field strength 
Bz,oo is characterized by the dimensionless flux-to-mass 
ratio To = B^^oo/pTrG^/^So.rof]- The reference clouds are 
not in an equilibrium state and are allowed to evolve into 
one with magnetic field £rozen-in, before ambipolar dif- 



fusion is turned on at time t = Q. At i = 0, we im- 
pose on top of the equilibrium surface density So(a;,2/) an 
m = 2 perturbation of relative amplitude A, T,{x,y) = 
I]o(a;,y) [1 -I- A cos(2(/))], where cj) is the azimuthal angle 
measured from the a:;-axis. A slow rotation is also added 
to the cloud at this time, according to the prescription 
Tr / N 4 w r 

y<t>{x,y) = ==Cms, (2) 

^0 + V ^0 + 

where c„is = c,(l + Fq)^/^ is essentially the magnetosonic 
speed, and lo measures the rotation rate. The subse- 
quent, ambipolar diffusion-driven evolution of the per- 
turbed cloud is followed numerically, subject to the condi- 
tion of fixed p and B^ at the cloud outer radius, taken to 
be twice the characteristic radius tq. Other choices of the 
reference state will be explored elsewhere. 

Our calculations are carried out using an MHD code 
based on that of Nakamura & Hanawa (1997). The hydro- 
dynamic part is solved using Roe's method. The equation 
for magnetic field evolution has a form identical to that 
for mass conservation, and is solved in the same man- 
ner. Both the magnetic and gravitational potentials sat- 
isfy the Poisson equation, and are solved using a convolu- 
tion method (Hockney & Eastwood 1981). We refine the 
computational grids whenever the so-called Jeans condi- 
tion is about to be violated (Truelove et al. 1997), tak- 
ing into account the magnetic effect on Jeans condition. 
The computations are carried out using nondimensional 
quantities. The units we adopted are Cg for speed, So,rcf 
for surface density, 27rGEo.rof for gravitational accelera- 
tion, and i?2,oo for magnetic field strength. The units for 
length and time are, respectively, Lq = Cg/(27rGEo,ref) and 
to = c,/(27rGEo ,ref / * 

3. NUMERICAL RESULTS 

The results of our calculations are shown in two fig- 
ures. In Figure 1, we show the evolution of a "benchmark" 
cloud with the dimensionless radius ro = 7.57r, initial fiux- 
to-mass ratio Fq = 1.5 and rotation rate ix> = 0.1. The 
m = 2 perturbation added has a relative amplitude of 
"merely" A = 0.05. During the initial, quasi-static stage 
of evolution, a central core condenses gradually out of the 
magnetically subcritical cloud, with no apparent tendency 
for the mode to grow. Rather, the iso-density contours 
appear to oscillate, changing the direction of (slight) elon- 
gation along the .T-axis in the disk plane to y-axis, as ev- 
ident from panels (a) and (b). (We confirmed that the 
disk continues to oscillate around the equilibrium state in 
the absence of ambipolar diffusion.) The period of the 
oscillation is more or less independent of a; as long as to 
is small. After a supercritical core develops, the contrac- 
tion becomes dynamic, as seen from the velocity field in 
panel (c). By this time, the m = 2 mode has grown sig- 
nificantly, resulting in a bar-like core at the center with 
an aspect ratio of roughly 2. As the collapse continues, 
the bar growth remains relatively modest [panel (d)] until 
the very end of the starless collapse, when the growth rate 
increases dramatically [panel (e)]. Beyond a critical di- 
mensionless surface density of 1.9 x 10^, corresponding to 
a volume density of nu — 10^^ cm^^, we change the equa- 
tion of state from isothermal to adiabatic, with an index 
of 5/3, to mimic the transition to optically thick regime 
of cloud evolution. The density distribution and velocity 
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field around the time when an accretion shock begins to 
develop are shown in panel (f). The bar bound by the 
shock is analogous to the "first" core of spherical calcu- 
lations. It has an aspect ratio of about 13. This highly 
elongated "first" bar is expected to break up into two or 
more pieces, as discussed further in § 4. 

To discuss the bar growth more quantitatively, we plot 
in Figure 2 the aspect ratio R of the darkest region in 
Figure 1 (with S > 10~^/^Sniax) as a function of the cen- 
tral surface density Smax: together with the flux-to-mass 
ratio r at the cloud center. It is clear that the m = 2 
mode grows little when the entire cloud remains subcrit- 
ical (i.e., F > 1), as one might expect intuitively. As the 
central region becomes more and more supercritical, the 
mode grows significantly, producing a bar of aspect ra- 
tio ~ 2. This ratio stays more or less "frozen" during 
the subsequent evolution, presumably due to the onset of 
rapid collapse, which leaves little time for the bar to grow. 
The flux-to-mass ratio does not change much either dur- 
ing this period, also as a result of the rapid collapse which 
prevents the magnetic flux from leaking out. The fairly 
strong trapped field (with F ^ 0.5) further hinders the 
bar growth, which picks up speed only towards the end of 
collapse before forming a singularity, as a result of the Lin, 
Mestel, & Shu (1965) instability. 

We have studied other models with different initial pa- 
rameters in the range of 1.25 < F < 2, 57r < ro < IOtt, 
< uj < 0.1, and 0.05 < A < 0.1. We find that the bar 
grows more rapidly for a smaller Fq, larger ro, larger w, or 
larger A. For example, if we choose Fq = 1.25, ro = 7.57r 
(adopted by Ciolek & Basu (2000) for their axisymmetric 
model of the well-observed starless core L1544), ui = 0.1, 
and A = 0.05, then the aspect ratio reaches 33 when the 
"first" bar forms, as shown in Figure 2. Interestingly, a 
plateau region is also evident, where the aspect ratio is 
temporarily "frozen" at i? 2, as in the benchmark model. 
We find that the frozen value of the aspect ratio depends 
weakly on the initial amplitude of perturbation, which may 
have observable consequences as we discuss below. 

4. DISCUSSION AND CONCLUSION 

We have followed numerically the growth of to = 2 mode 
during the ambipolar-diffusion driven evolution of mag- 
netically subcritical clouds. Our main conclusion is that, 
despite the presence of the strong magnetic field, a per- 
turbation of modest amplitude can grow nonlinearly into 
a highly elongated bar, which is expected to fragment into 
small pieces gravitationally. Before discussing fragmen- 
tation, we comment on the effects of bar growth on the 
observed shapes of molecular cloud cores. 

Dense cores of molecular clouds are intimately associ- 
ated with star formation, with roughly half of them al- 
ready harboring infrared sources (Beichman et al. 1986). 
The other half are thought to be well on their way to star 
formation, with the "starless" phase lasting for only a few 
dynamic times. The cores are observed to have significant 
elongation, with a typical aspect ratio of 2 (Myers et al. 
1991). Statistics of core elongation have been interpreted 
as indicating most cores have prolate 3-dimensional shape 
(Ryden 1996; Curry & Stabler 2001), although oblate cores 



formed as a result of settling along magnetic field lines can 
have a projected aspect ratio similar to that observed (Li 
& Shu 1996). Nonlinear growth of the bar mode in a di- 
rection perpendicular to the field lines modifies the core 
shape drastically, making it triaxial in general. Some evi- 
dence for the triaxial nature of cores has been marshalled 
by Basu (2000). We propose that it is due to the significant 
bar growth during the transition period when F decreases 
from ~ 1 to ~ 0.5, after the core has become supercritical 
(which makes bar growth possible) but before a rapid col- 
lapse sets in (which leaves little time for the bar to grow). 
The bar growth removes an oft-invoked objection against 
the standard scenario of ambipolar diffusion-driven core 
formation based on the elongated or filamentary shape of 
dense cores (Ward-Thompson, Motte & Andre 1999). In- 
deed, the triaxial nature of the dense cores, coupled with 
the relatively slow, subsonic infall motion inferred (Lee, 
Myers & Taffala 2001), may provide the strongest support 
yet for the standard scenario. 

The continued evolution of bar-like cores, such as the 
one shown in Figure Ic, would formally lead to a singular 
filament, if the isothermal assumption is kept. However, 
it is well known that the equation of state stiffens when 
the optical depth exceeds unity (see however Masunaga 
& Inutsuka 1999). The stiflfening slows down the collapse, 
allowing more time for the elongated bar to fragment grav- 
itationally into small pieces. In fact, the major axis of the 
bar exceeds twice the critical wavelength for fragmenta- 
tion of an infinitely- long filament (Larson 1985) in pan- 
els (e) and (f) of Figure 1. Such bar formation is also 
shown by Matsumoto & Hanawa (1999) and Sigalotti & 
Klapp (2001) for nonmagnetized clouds. Recent examples 
of (nonmagnetic) bar fragmentation induced by the stiffen- 
ing of equation of state are given in Boss et al. (2000) and 
Sigalotti & Klapp (2001). For magnetized bars, one needs 
to consider in addition the rapid decoupling of magnetic 
field from matter, which occurs above a density of order 
10^° cm-3 (Nishi, Nakano, & Umebayashi 1991). The de- 
coupling decreases the magnetic support quickly, and is in 
some sense equivalent to a sudden cooling. It should make 
the fragmentation easier. 

Boss (2000) studied the fragmentation of 3D magnetic 
clouds numerically, treating the magnetic forces and am- 
bipolar diffusion in an approximate way. He concluded 
that magnetic fields can enhance cloud fragmentation, by 
reducing the tendency for the development of a central sin- 
gularity, which would make fragmentation more difficult. 
We also find that magnetic fields can promote fragmenta- 
tion, but for a different reason. Strong magnetic fields can 
support clouds of more than one Jeans mass, which pro- 
vides an initial condition that is conducive to fragmenta- 
tion once the magnetic support weakens. We believe that 
it is the multi-Jeans mass nature of the magnetically sub- 
critical clouds that drives the nonlinear growth of the bar 
(and higher order) mode. We conclude that the magneti- 
cally subcritical clouds envisoned in the standard scenario 
of star formation can produce not only single stars but 
also, perhaps even preferentially, binary and multiple star 
systems. A parameter survey will be carried out to firm 
up this conclusion. 
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Fig. 1. — Surface density and velocity distributions for the model with (ro,ro,ti-') = (1.5, 7.57r, 0.1) at six dimensionless times; (a) t=0, 
(b) 17.49, (c) 24.46, (d) 25.21, (e) 25.2500, and (f) 25.2503. In panels (c) through (f), only the central regions are shown. The contours 
in each panel are for surface density normalized by the maximum, with values indicated by the grayscale. The numbers above each panel 
are the maximum surface density (Smax), fiux-to-iiiass ratio (F) at the density peak, and length unit for each panel. The arrows show the 
velocity vectors and are normalized by the arrow above each panel. If we choose Sq ref = 0.01 g cm~^, Tgfj = 30 K, then the units for time, 
length and speed would be to = 2.46 x 10^ yr, Lo = 0.082 pc and Cs = 0.33 km s~^. The dimensionless flux-to-mass ratios averaged in the 
darkest regions of the panels are (a) 4.46, (b) 2.77, (c) 1.26, (d) 0.645, (c) 0.503, and (f) 0.480, respectively. Note that the "bar" as defined 
by the darkest region becomes completely supercritical at late times, even though the cloud as a whole remains subcritical. The maximum 
volume densities in panels (a) through (f) would be (a) 4.46 x 10^ cm-3, (b) 1.70 x lO'' cm-3, (c) 2.79 x 10^ cm-3 (d) 2.80 x lO'^ cm-3, (e) 
2.78 X lO^i cm-3, and (f) 4.77 x lO^^ cm-3, respectively. 



Bar Formation in Magnetic Clouds 



5 




Fig. 2. — Evolution of the aspect ratio R of the bar (solid lines) and the flux-to-mass ratio F (in units of the critical value) at the surface 
density maximum (dashed lines). Thick and thin lines are for models with (Fo,ro,t<^) = (1.50, 7.5tt, 0.1) and (ro,ro,ijj) = (1.25, 7.57r, 0.1), 
respectively. For reference, lines of F = 1 and R = 1 are also plotted. 



